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■^r ' We develop a combined hydro-kinetic approach which incorporates a hydrodynamical ex- 

pansion of the systems formed in A+A coUisions and their dynamical decoupling described by 

^i^j ' escape probabilities. The method corresponds to a generalized relaxation time (tj-ci) approxi- 

mation for the Boltzmann equation applied to inhomogeneous expanding systems; at small t^cI 
it also allows one to catch the viscous effects in hadronic component - hadron-rcsonance gas. 
We demonstrate how the approximation of sudden freeze-out can be obtained within this dy- 
namical picture of continuous emission and find that hypersurfaces, corresponding to a sharp 
y I freeze-out limit, are momentum dependent. The pion nix spectra are computed in the devel- 

oped hydro-kinetic model, and compared with those obtained from ideal hydrodynamics with 
the Cooper-Frye isothermal prescription. Our results indicate that there does not exist a univer- 
sal freeze-out temperature for pious with different momenta, and support an earlier decoupling 
of higher px particles. By performing numerical simulations for various initial conditions and 
equations of state we identify several characteristic features of the bulk QCD matter evolution 



f^ , preferred in view of the current analysis of heavy ion collisions at RHIC energies. 
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1 Introduction 

Hydrodynamic models successfully describe basic features of high energy nuclear collisions at CERN 
SPS and especially at BNL RHIC (for reviews see, e.g., Refs. [US [3]), where the utilization of 
ideal hydrodynamics was supported by theoretical results: it was advocated [4] that deconfined 
matter behaves like a perfect liquid. The hydrodynamic approach to A-\-A collisions implies a 
very hot and dense strongly interacting matter as the initial state. Such a state is assumed to be 
formed soon after the collision, then expands hydrodynamically until the stage when the picture 
of the continuous medium is destroyed. Roughly, it happens when the mean free path of particles 
becomes comparable with the smallest characteristic dimension of the system: its geometrical size 
or hydrodynamic length. This approach allows one to account for the complicated evolution of 
the system at a preconfined stage and in the vicinity of the possible phase transitions by means 
of a corresponding equation of state (EoS). By studying the spectra of different particle species 
versus the initial conditions and EoS, one could get information on the thermal partonic stage of 
the "Little Bang" and further system evolution, which could be also conclusive in the context of 
discriminating the possible phase transitions. 



The problem is, however, whether the predicted momentum spectra, obtained in hydrodynamic 
models for given initial conditions and EoS, are unambiguous. These spectra depend not only 
on the initial conditions but also on the final conditions of hydrodynamic expansion since one 
cannot use hydrodynamic equations for infinitely large times because the resulting very small 
densities destroy the picture of continuous medium of real particles. Evidently, the simplest receipt 
of the spectrum calculation is the Cooper-Frye prescription (CFp) [5] which ignores the post- 
hydrodynamic (kinetic) stage of matter evolution and assumes that perfect fluid hydrodynamics is 
valid till some 3D hypersurface, e.g., as was supposed by Landau [6], till the isotherm T ~ m.,^, 
where sudden transition from local thermal equilibrium to free streaming is assumed. 

It has been known for a long time that CFp leads to inconsistencies [Zl [8], if the freeze-out 
hypersurface contains the non-space- like sectors, and should be modified to exclude formally nega- 
tive contributions to the particle number at the corresponding momenta. The simplest prescription 
is to present the distribution function as a product of a local thermal distribution and the step 
function like 9{pfj,n^{x)) [8], n^Ufj, = ±1, where n^ is a time-like or space-like outward normal to 
a freeze-out hypersurface a. Thereby freeze-out is restricted to those particles for which Pfj_n^{x) 
is positive. This receipt was used in Ref. [9] to describe particle emission from enclosed freeze-out 
hypersurface with non-space-like sectors and it was found that a satisfactory description of central 
Au+Au collisions at RHIC was reached for a physically reasonable set of parameters. The main 
features of the experimental data were reproduced: in particular, the obtained ratio of the out- 
ward to sideward interferometry radii is less than unity and decreases with increasing transverse 
momenta of pion pairs. Thereby, the results of Ref. [9] clearly indicate that early particle emission 
off the surface of the hydro dynamically expanding fireball could be essential for proper description 
of matter evolution in A+A collisions. 

However, sharp freeze-out at some 3D hypersurfaces is a rather rough approximation of the 
spectrum formation, because the particle emission process of fireballs created in high energy heavy 
ion collisions is gradual in time. Results of many studies based on cascade models contradict the idea 
of sudden freeze-out and demonstrate that in fact particles are emitted from the 4D volume during 
the whole period of the system evolution, and deviations from local equilibrium conditioned by 
continuous emission should take place (see, e.g., [TO]). Moreover, freeze-out hypersurfaces typically 
contain non-space-like parts that lead to a problem with energy-momentum conservation law in 
realistic dynamical models [8]. This concerns also hybrid models [H] where the transport model 
matches hydrodynamics on such a kind of isothermal hypersurface of hadronization [12] |^ 

An attempt to introduce 4D continuous emission in hydrodynamic framework has been done in 
Ref. [13j within a simple ideal hydrodynamics, and also supported by numerical transport codes 
calculations, e.g., Ref. [H]. It was found in these papers that the transport freeze-out process 
is similar to evaporation: high-pT particles freeze out early from the surface, while low-pT ones 
decouple later from the system's center. 

The idea of continuous decoupling was further developed based on Boltzmann equations in 
Ref. [15j where, in particular, an approximate method that accounts for the back reaction of 
the emission on the fluid dynamics was proposed. It is worth noting that the back reaction is 
not reduced just to an energy- momentum recoiling of emitted particles on the expanding thermal 



^ Note that in this hybrid picture the initial conditions for hadronic cascade calculations could be formulated 
also on some (arbitrary) space-like hypersurface where, however, hadronic distributions are deviated from the local 
equilibrium, in particular, because of an opacity effect for hadrons which are created during a "mixed" stage of phase 
transition. These nonequilibrium effects could seriously influence the results of hybrid models in its modern form 



medium, but also leads to a rearrangement of the medium, producing a deviation of its state from 
the local equilibrium, accompanied by changing of the local temperature, densities and collective 
velocity field. This complex effect is mainly a consequence of the impossibility to split the evolution 
of the single finite system of hadrons into the two components: expansion of the interacting locally 
equilibrated medium and free streaming of emitted particles, which the system consists of. Such a 
splitting, accounting only for the momentum-energy conservation law, contradicts the underlying 
dynamical equations such as a Boltzmann one [l5j. In view of this, the ideas proposed and results 
obtained in the quasiclassical approach should be a clue for a quantum treatment of the problem 
which recently begins to to be developed [16] o Note, however, that the estimates of the influence of 
the quantum effects, such as the distortion of the wave function, on the spectra and Bose-Einstein 
correlations testify to be relatively small corrections to the quasiclassical approximation [18j . 

The aim of the present work is twofold. First, we develop the formalism of the hydro-kinetic 
approach and propose approximations for practical calculations. The problem of self-consistency of 
the method accounting simultaneously for both the particle emission and fluid (in general, viscous) 
dynamics is studied in detail. Special attention is also paid to the discussion of the applicability 
conditions of the Cooper- Frye prescription for sharp freeze-out. 

Second, we develop a simple hydro-kinetic model of hadronic emission that describes the evo- 
lution and emission of particles from a hydrodynamically expanding system undergoing a phase 
transition. For the sake of simplicity we consider here one type of escaping particle (pions) only 
and do not take into account in numerical calculations the back reaction of the emission on fluid 
dynamics. We study within such an approach what type of initial conditions, equation of state, 
etc., are preferred in view of current analysis of heavy ion collisions at RHIC energies. 

2 Hydro-kinetic formalism for heavy ion collisions 



It was proposed in Ref. [15] to describe the hadronic momentum spectra in A+A collisions based 
on the escape function of particles which are gradually liberated from hydrodynamically expanding 
systems. The escape function, introduced in [13], is calculated within the Boltzmann equations in 
a specific approximation based on a hydrodynamic approach. It was shown that such a picture 
corresponds to a relativistic kinetic equation with the relaxation time approximation for the collision 
term, where the relaxation time tends to infinity, r^ei — > oo, when t ^ oo, indicating a transition 
to the free streaming regime. For one component system the equation has the form: 

pf" df{x,p) ^ f{x,p) -/^q(a;,p) 
Po ^x^' Treiix,p) 

Here f{x,p) is the phase-space distribution function, /^ ^^{x,p) is the local equilibrium distribution 
with local velocities, temperatures, and chemical potentials that should be found from Eq. ([T]) and 
the initial fo{x,p), and Tj-(,i{x,p) is the relaxation time (inverse rate of collisions in gases), 

rUx,p) = '-^^^. (2) 

Here r*|,[(x,p) is related to the local fluid rest frame (local rest frame of the energy flow) where the 
collective four-velocity is u^ = (1,0). The relaxation time depends on the cross section and is a 
functional of /' ^^{x,p). 



^See also the very recent contribution to this field in Ref. [17) . 



As it is well known |191 I20j . such kinds of equations at Trei <C Texp (inverse of expansion 
rate) describe in the first approximation the viscosity effects in gases with a coefficient of shear 
viscosity r] oc Trel?^^. Therefore the method explained below catches in the first approximation also 
the viscosity effects in an expanding hadronic gas, characterized by fields of temperatures T and 
particle densities n. The viscosity effects in the quark-gluon plasma (QGP) evolution cannot be 
described in this way because strongly interacting QGP is not a gas, but almost an ideal liquid [3]. 

The formal solution of Eq. ([T]) can be presented in the following form: 

f(t,r,p) = f(to,r {t-to),p)exp\- / p— — -ds > + 



, PO 

to 



rrei(t',r-^(t-t'),p) """^^ 
to 



T-rel{s,r-^{t-s),p) 



ds 



> dt', (3) 



where f{to,r,p) is the initial distribution at t = to- The relaxation time r*g[ as well as the local 
equilibrium distribution function /' '^^ are functionals of hydrodynamic variables: temperature T, 
chemical potential //, and collective four- velocity n^. The space-time dependence of the corre- 
sponding variables are determined by demanding the local conservation of the energy-momentum 
with tensor T^'^{x) and, if necessary, net particle number, with current n^{x) (assuming no particle 
production) 

d^T^^'^ix) = 0, (4) 

d^n^ix) = 0, (5) 

where (see, e.g., Ref. [20]) 

—k^'^f{x,k), (6) 

/d^k 
—k^f{x,k). (7) 



These conservation laws lead to rather complicated equations for hydrodynamic variables. It is 
worth noting that for an expanding system the relaxation time T*^i{x,p) increases with time and, 
therefore, the deviations from local equilibrium increase too, thereby preventing a use of the widely 
applied approximate methods based on the expansion of the distribution function in the vicinity of 
the local equilibrium. 

Then to solve the kinetic equation ([T]), in accordance with the conservation laws ([1| and ^, 
we need an approximate method that could be applied even for strong deviations from local equi- 
librium. It is not our aim here to suggest an exclusive solution of the problem. Rather some 
arguments are presented below by the example of the relativistic one component Boltzmann gas 
with particle number conservation, 

/^''(x,p) = (2vr)-3exp(-^^^^), (8) 

to show that such a method could be developed based on the following procedure. 



To take into account nonequilibrium effects accompanying the particle emission in inhomoge- 
neous violently expanding systems, we utilize the integral representation ^ of kinetic equation 
([1]). Then, performing a partial integration of the second term in Eq. ^ and, assuming that 
f{tQ,r,p) = f^^^{tQ,r,p), one can decompose the distribution function to a local equilibrium part, 
/ ^^, and a part describing a deviation from the local equilibrium behavior, g: 

f = f'^{x,p)+g{x,p), (9) 

where 

^ ""^ ' ' ' • ' }dt'. (10) 



g{x,p)=- / — exp<(- / p— — -ds 

J dt' \ J r,ci(s,r-^(t-s),p) 

to ^ t' '^° 

Note that both functions, / "^^ and g, are functionals of hydrodynamic variables, g depends 
also on the relaxation time r^ei that defines the mean time interval between collisions, and Tj-ei 
depends in its turn on the distribution function / ^^ and the cross section. The evolution of 
the distribution function f{x,p) should satisfy the energy-momentum conservation and, because 
r^^[/] = r^!^[/ieq + g] = j^M^fjieq] _^ T'"'[g] for systcms where the interaction energy can be 
neglected, it takes the form of hydrodynamic equations for the perfect fluid with "source" , 

d,T''''[f'^^]=G''[g], (11) 

where 

G%] = -d,T'^%]. (12) 

The equation that takes into account the conservation of particle number has a similar form: 

d,n'^[f''i] = S[g], (13) 

where 

S[g] = -d,n-[g]. (14) 

To find an approximate solution of Eqs. dHJ-dH]), one can solve the equations 

5^r^^[/'^^] = 0, (15) 

5.n''[/<=^] = 0, (16) 

and, thereby, utilize the hydrodynamic variables in the perfect fluid approximation. Namely, the 
hydrodynamic variables in this approximation can be used to calculate the deviation from local 
equilibrium g{x,p) according to Eq. (jlOp and, then, "source" terms G'^[fl'] and S[g\ on the right- 
hand sides of Eqs. (jlip and (|13p . Then the left-hand sides of these equations are functionals of local 
equilibrium functions and have the simple ideal fluid forms, while the right-hand sides associated 
with a "source" are explicit functions which describe a deviation from the local equilibrium and 
depend on hydrodynamic variables in the perfect fluid approximation: 

d^T'^^if ^^(T, u,, m)] = G^[Tid, <, Wd, r;fi], (17) 

d,n^[f ^^(T, n^, ^)] = S[rid, njf , ;Uid, r^^], (18) 



where, for one component Boltzniann gas with elastic colhsions only, the relaxation time r^l^j is 
the inverse of collision rate in ideal fluid, R {x,p), and has the following form (in the co-moving 
frame) : 



^ (^'P) = / T^^TTl^exp — — a{s)^^—-— . (19) 



iiii^^P) J (2^)' V Tid{x) J ' ' 2EpEk 



Here Ep = y^p^ + m?^ E^ = Vk^ + m-'^, s = {p + k)"^ is the squared cm. energy of the pair, 
and a{s) is the corresponding cross section. A solution (T(x),m^(j;),/x(x)) of Eqs. (fT7|l . (fTSll 
accounts for the back reaction of the emission process on hydro-evolution and provides us with the 
hydrodynamic parameters which finally should be used to calculate the locally equilibrated part 
/'°'^(x,p) of the complete distribution function f{x,p). Then, the distribution function obtained 
in this way, f{x,p) = /^ '^^[T, n^,/x] -|- g[Ti(j,ii)^,/iid,Tj!gj], satisfies the conservation laws, takes into 
account the nonequilibrium peculiarities of the evolution and is constructed in agreement with the 
corresponding EoS. Of course, this scheme allows us to make the next iterations in solving Eq. ([1]). 
Note also that because the "source" term on the right-hand side of Eq. (I17p is a known function, 
the causality is preserved in this description of dissipative systems. 

3 Kinetics of the freeze-out in hydro-kinetic approach 

The approach developed in the previous section allows us to study an important problem: under 
which conditions the Landau/Cooper-Prye prescription (CFp) of sudden freeze-out is valid and 
how, then, to define the corresponding 3D freeze-out hypersurface. The CFp is traditionally uti- 
lized when hydrodynamics is applied to describe the later stage of the matter evolution in A+A 
collisions, and also is a basic ingredient of various hydro- motivated parametrizations (see, e.g., 
[HlElj). A widely used heuristic freeze-out criterion for finding the freeze-out hypersurface of a 
violently expanding system is either the comparability of the hydrodynamic rate of expansion with 
the kinetic rate of collisions [22] or of the mean free path of particles with the geometrical size of 
the system [6l[23]. While the problem was extensively studied before (see, e.g., Refs. 13 O [151 [23]), 
a complete understanding of the conditions allowing the utilization of the sharp freeze-out approx- 
imation (that corresponds formally to sudden transition from local equilibrium to free streaming) 
of spectra formation and an unambiguous definition of the corresponding freeze-out hypersurface 
are still absent. A discussion in this section based on the analytic approximation to a calculation of 
momentum spectra could be, in our opinion, useful for understanding the conditions of applicability 
of the CFp and improvement of it, if necessary. 

Let us integrate over space variables distribution ([3|) to represent the particle momentum density 
at large enough time, t — > oo, when particles in the system stop to interact: 

j3 AT f 

^{t)^n{t,p) = J d^rf{t,r,p). (20) 

The result can be presented in the general form found in Ref. [15j : 

n{t -^ oo,p) = / (f'rf{to,r,p)'P{to,r,p) + 
t 
I d^r f dt'F^'''''{t',r,p)V{t',r,p). (21) 

to 



The term F^^^'^{x,p) corresponds to "gain" term in the Boltzmann equation and is associated with 
the number of particles which came to the phase-space point {x,p) just after the interaction in the 
system. The probabihty V{x,p) for a particle to escape the system from space-time point x = (t', r) 
is expressed explicitly in terms of the rate of collisions R along the world line of the free particle 
with momentum p, 

V{t',r,p)=exp\- [dsR{s,r + ^{s-t'),p)] . (22) 

V / ^° J 

The first term in Eq. (|2ip describes the contribution to the momentum spectrum from particles that 
are emitted from the very initial time, while the second one describes the continuous emission with 
emission density S = F^^^'^V from 4D volume delimited by the initial and final (where particles 
stop to interact) 3D hypersurfaces [15]. In a particular case, which we consider in this article, 
R{x,p) = l/rrei(x,p) and FS''™(x,p) = f^''\x,p)/T^e\{x,p). 

To provide straightforward calculations leading to the Cooper-Frye approximation let us shift 
the spacial variables, r' = r + ^(to — t')? ™ the second term of Eq. (j2ip aiming to eliminate the 
variable t' in the argument of the function R = l/r^ei in (|22p . Then 



n{t,p) = / cfrf{to,r,p)exp < - / : — r-t — ^ — -ds \ + 



1 

T-Tciis,r + ^{s-to),p)' 



, Po 

to 



fd^r' f dt'f''i{t',r' + ^{t' -to),p)Q{t',r',p), (23) 

J J Po 

to 
where 

Here, / '^^, Trei are functionals of temperature, hydrodynamic four-velocity, and chemical potential, 
and these quantities are governed by hydrodynamic equations (|lip -(|14p. Note that 

Q{t',r',p) = -^P{t',r',p), (25) 

where P{t',r',p) is connected with the escape probability V(t',r,p): 

P{t',r',p) = exp}- f 



t 

1 

-ds 



T,ci{s,r' + f-{s-to),p) 
t' 



^=p(t',r'+ P-(t'-to),p). (26) 

Po 



Therefore 

oo 

' dt'Q{t',r',p) = l-V{toy,p). (27) 



/ 



to 
The 4D emission density has the form 



S{t',r',p) = f'^{t',r' + ^{t' -to),p)Q{t'y,p). (28) 

Po 
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In order to have a tractable approach and set up the conditions of vahdity of CFp, let us assume 
that at each r' and p there is a maximum in t' of the emission function S{t' , r' ,p) inside the interval 
[to,t\, and that the position of the maximum, t' = t'^{r',p), is mainly conditioned by Q{t',r',p). 
Then corresponding hypersurface t'^{r',p) is defined by the conditions 

dQ{t',r',p) _ 

— dF' — 1*'=*^-°' (29) 

^2 \t'=t'. < 0' (30) 

and, utilizing Eq. ([M|) . we get 



dT,,i{t',r' + ^^{t'-to),p)\t'=t'^ 

dt' 
d^TUt',r' + ^{t'-to),p)\t'=t'^ 



= 1, (31) 

> 0, (32) 



where d/dt' is the full time derivative. 

Let us demonstrate that Eqs. (j3ip . (j32p generalize the heuristic freeze-out criterion [22]. Ac- 
cording to the latter, the freeze-out happens when 

Tscat ~ '''cxp ) V""^/ 

where Tgcatix) is the mean time interval between successive scattering events, 

rscat(x) = - — . . . . , (34) 

{va)[x)n{x) 

V is the relative velocity between the scattering particles, n(x) is the particle number density and 
o" is the corresponding total cross section and the sharp brackets mean an average over the local 
thermal distribution. The inverse hydrodynamic expansion rate Texp(2;) is the collective expansion 
time scale, 

--<-' = -(i) ""«-")" = -(;i)l^)"> P^' 

where t* is the proper time in the local rest frame of the fluid. 

To demonstrate that the heuristic freeze-out criterion (f33l) follows from Eqs. pT]) and (f32]) 
under certain conditions, we note first that Tsca.t{x) is similar to t^^*, see Eq. (J19p . Then, rewriting 
Eq. (I3ip in the form 

■^rel(t',r,p)|t/=tj^^r=r'+J^(t^-to) = 1' (36) 



d p d 
dt' pq dr 

we get from the above equation that 



1 \~V 5 P d\ / 1 



+ — — I I r I \t'=t'^,r=r'+^{t',-to) 



r,ei(t^,r,p); \dt' po 9r 7 V^rei(t',r,p)y "-*-"-'■ +TO* 

1 

r.el(t'.,r,p)"*'=*'^''--'+^(*-*o) 



(37) 



Now, if we turn to the local rest frame of fluid and neglect there the momentum dependence of 
the above equation, as well as deviations from local equilibrium, we get, accounting for Eqs. ()19p 
and (pll) . that the left-hand side of Eq. (f37|l is approximately — (^^) \t'*=t'* while the right- 
hand side is 1/T^^*{t'*,r*) w 1/Tsca.t{t'^,''^*), recovering thereby the criterion ([33]) of the freeze-out. 
One can note, however, that unlike the heuristic definition ()33p the true freeze-out hypersurface 
t'^{r',p) depends on momentuno and thereby particles of different momenta freeze out on different 
hypersurfaces. This result has already been observed in [13] in terms of escape probabilities. 

Now, we proceed to the Cooper-Frye representation of the freeze-out process. For this aim it is 
convenient to introduce a new variable in the second term of Eq. (I23p . namely, 

(38) 



(39) 



r = r 


PC 


-{t'^{r',p)-to). 

) 


Then 


r' = r {ta{r,p) -to), 

Po 


and the corresponding Jacobian is 


det 




_ . P dt„ 
Po dr ■ 



(40) 

This function is positive as will be shown below. Then Eq. (j23p takes the form 

n(i -^ oo,p) = / (frf{to,r,p)'P{to,r,p) + 
t 

I d'r (l - P-^) / dt' f-\t\r + ^{t' -t,),p)Q{t\v - ^{t,{r,p) -to),p). (41) 

J \ Po or J J Po Po 

to 

To estimate the accuracy of the Cooper-Prye approximation let us apply the saddle-point method 
for the calculation of the second term in Eq. ()4ip by the use of expansion of the function lnQ{t' ,r — 

£-{t^{r,p)-to),p) =Hr-l{t',r + ^{t' -t,{r,p),p)V{t',r+^{t' -t^{r,p),p)] (see Eq. dMD) near 
the point of maximum in t': 

oo 

exp(- / (is;-(^^^rp^(7i:7;(7;^) 

Q{t',r-^{U{r,p)-to),p)^ '- .. , ° , exp{-it' -t^{r,p)f/2D\r,p)),{42) 

Po Tj.ei{t^{r,p),r,p) 

where 

n2/ N r,ci(t^(r,p),r,p) 

dt^ \t'=U{r,p),r'=v-^{U{v,p)-to) 

Let us suppose that probability V{to,r,p) for particles to escape just at the initial moment is 
negligible, except, of course, the periphery of the system at t = to which, however, gives a relatively 



The momentuin dependence of freeze-out hypersurface is conditioned by both momentum and spacial dependen- 
cies of Trei- Note that the momentum dependence of the relaxation time was explicitly demonstrated recently in Ref. 
[25) based on the numerical results of parton cascade simulations. 



small contribution to the spectra. Then, taking into account that in the region of freeze-out the 

oo 

particle suffers the last collision and, thus, / dsr^^ (s,r + ^(s — ta{r,p),p) ~ 1 {t~^ is the rate 

of collisions), and also accounting for the normalizing condition (j27p for Q, one can integrate Eq. 
over t' and estimate the temporal width of the emission zone 



Dir,p) = —=Trei{ta,r,p)exp{ / dsT^J{s,r + —{s - tair,p),p) \ ~ Trci(ta, r,p)- 



(44) 



The above relation means that at the freeze-out hypersurface Treir^'gi — -*■' ^^^ ^^• 

Therefore if the temporal homogeneity length A(r,p) of the distribution function / °^ near the 
four-point {ta{r,p),r) is much larger then the width of the emission zone, A(r,p) >> Tj-c\{r,p), then 
one can approximate /^ ^'^{t', r+-^{t' — tf^),p) by /^ '^^{ta, r, p) in Eq. (j4ip and perform integration 
over t' accounting for normalizing condition (j27p . As a result we get momentum spectra in the 
Cooper-Frye form: 

p\it -^oo,p) = J da^p^'f^^ix^p), (45) 

where 

da^p^ = d^r L - P^) . (46) 

Now let us summarize the conditions for the utilization of the Landau/Cooper-Prye approxi- 
mation of sudden freeze-out. They are the following: 

(i) For each momentum p, there is a region of r where the emission function as well as the 
function Q (see Eq. (I24p ) have a sharp maximum with temporal width D(r,p). The formal 
conditions of the maximum are defined by Eqs. (|29p . (j30p . 

(ii) The width of the maximum, which in the case of one component system is just the relaxation 
time (inverse of collision rate), should be smaller than the corresponding temporal homogeneity 
length of the distribution function: A(r,p) » D{r,p) ~ ri.ci(r,p). 

(iii) The contribution to the spectra from the residual region of r where the saddle point method 
(Gaussian approximation (j42p and/or condition Tj-d <^ A) is violated does not affect essentially the 
particle momentum density. In other words, the space regions where for fixed momentum p either 
there is no sharp maximum of the emission function or phase-space density changes very rapidly 
are not important for spectra formation. Note, that this condition is most questionable and has to 
be checked in realistic dynamical approaches for freeze-out, like hydro-kinetics, or transport models 
with appropriate initial conditions for hadronic phase. 

(iv) The escape probabilities V{tQ,r,p) for particles to be liberated just from the initial hyper- 
surface io are small almost in the whole spacial region (except maybe peripheral points) and so, 
one can neglect the first integral in Eq. ()4ip . 

If the conditions (i) - (iv) are satisfied, then the momentum spectra can be presented in Cooper- 
Frye form in spite of the fact that there is no sudden freeze-out and the decaying region has a 
finite temporal width D{r,p) ~ Trei{t^(r,p),r,p). Also, what is very important, such a generalized 
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Cooper- Frye representation is related to the freeze-out hypersurface that depends on momentum p 
and does not necessarily enclose the initially dense matter (it will be demonstrated in the next 
section by numerical calculations of emission functions). 

In as much as the hypersurface to-(r,p) corresponds to the maximum of the emission of the 
particles with momentum p into vacuum, these particles cannot be emitted from point {ta{T,p),r) 
inward of the system. In the latter case neither the emission function S nor function Q can have 
a maximum at this phase-space point; in fact, their values are near zero. Even formally, in the 
Gaussian approximation (j42p for Q, validated in the region of its maximal value, the integral 

oo 

/ ^'^'^rci (^^^ + ^(^ ~ ^o-(r,p),p) ^ 1, if the particle world line crosses almost the whole system. 

The latter results in Q ^ and, therefore, completely destroys the saddle-point approximation 
()42p for Q. If the particle crosses some non-space-like part of the hypersurface a moving inward 
of the system, it corresponds to the condition p^da^ < [7j. Hence always the value p^da^{p) in 
the generalized Cooper-Frye formula (I45p is positive: p^da^{p) > across the hypersurface where 
a fairly sharp maximum of the emission of particles with momentum p is situated; and so the 
requirement p^da^{p) > is a necessary condition for t(j{T^p) to be a true hypersurface of the 
maximal emission. It means that hypersurfaces of maximal emission for a given momentum p may 
be open in space-time, not enclosing the high-density matter at initial time io) and different for 
different p. In other words, the momentum dependence of the freeze-out hypersurface, defined by 
Eqs. ()3ip and (j32l) . naturally restricts the freeze-out to those particles for which p'^d(Tfj_{r,p) is 
positive. 

Therefore, there are no negative contributions to the particle momentum density from non- 
space-like sectors of the freeze-out hypersurface, that is a well-known shortcoming of the Cooper- 
Frye prescription [71 18]; the negative contributions could appear only as a result of utilization of 
improper freeze-out hypersurface that roughly ignores its momentum dependence and so is common 
for all p. If, anyhow, such a common hypersurface will be used, e.g., as the hypersurface of the 
maximal particle number emission (integrated over p), there is no possibility to justify the approx- 
imate expression for momentum spectra similar to Eq. (j45p . Also, in that case there is no common 
phenomenological prescription, based on Heaviside step functions 6(pfj_da^), which allows one to 
eliminate the negative contributions to the momentum spectra when p'^da^ < 0. The prescription, 
proposed in [7], eliminates the negative contributions in the way which preserves the number of 
particles in the fluid element crossing the freeze-out hypersurface related to the maximum of total 
(averaged over p) particle emission. Therefore it takes into account that at the final stage the 
system is the only holder of emitted particles. Another prescription ^ ignores the particle number 
conservation considering decaying hadronic system rather as a star - practically unlimited reservoir 
of emitted photons/particles. Both prescriptions have a problem with momentum-energy conserva- 
tion laws at freeze-out. The approach developed here overcomes all above-mentioned problems by 
considering the continuous dynamical freeze-out that is consistent with Boltzmann equations and 
conservation laws. 

One can see from Eq. ()42p that in the general case with nonzero thickness of the emission 
layer, D ~ r^ei 7^ 0, the emission density S (j28p cannot be approximated by means of the local 
equilibrium distribution function /' ^^ smeared by a (proper) time factor exp(— (r — ro)^/5r^) with 
constant tq and 5t'^. Such an ansatz cannot be used in place of the proper emission function as it 
is often utilized in hydro- inspired parametrizations (see, e.g., Ref. [21]) to take into account the 
gradual character of the freeze-out process in heavy ion collisions. In fact, parameters tq and (5t^ 
should be dependent at least on the spacial coordinates, as it is explained in detail in Refs. [9|ll5j. 
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Equations (j42p and (|43p demonstrate that they depend on particle momentum as well and that the 
emission is not locally isotropic. 

4 Pionic emission in specific hydro-kinetic models 

In this section we present and discuss the results of numerical calculations obtained in a few concrete 
realizations of the hydro-kinetic model (HKM) . For the sake of simplicity we consider the emission 
of only one particle species (negative pions, vr") from the expanding fireball. Also we employ 
here the ideal fluid approximation for hydrodynamic variables in the integral representation ([3]) to 
calculate the emission densities and momentum spectra of pions. This approximation results in a 
little (less than 10% for our calculations) violation of the conservation laws. We describe, in all 
simulations, the locally equilibrated state of hadrons by relativistic Boltzmann distributions. 

The numerical results presented in this section were obtained on the basis of our original numeri- 
cal 3D ideal hydro-code that was developed based on the relativistic Godunov-type HLLE algorithm, 
described in detail in Ref. [26j. We use Bjorken coordinates ^27j r = i/t^ — z^ ^ rj = tanh~"'^(z/t), 
instead of Cartesian ones, the corresponding transformation of hydrodynamic equations and con- 
servative variables are described, e.g., in Refs. [21 [28]. A second order of accuracy in time of this 
algorithm is achieved by using the predictor-corrector scheme. To achieve second order of accu- 
racy in space we apply a linear distribution of conservative variables inside each fluid cell. Having 
determined the evolution in the ideal fluid approximation, we proceed in calculating the emission 
densities using numerical (tabulated) space-time dependencies of collective velocities and thermo- 
dynamic quantities. An 8- and 16-point Gaussian quadrature method implemented in ROOT [29| 
are used to calculate the escape probability for each space-time position and particle momentum. 
To calculate the resulting 4D integrals for particle spectra we use Monte Carlo method and check 
convergence of the results with increasing number of sampling points. We use ROOT for the results 
output also. 

First, we present our results on a specific toy-HKM that describes expanding one component 
relativistic ideal Boltzmann gas (p = nT) of vr^ with presumed conservation of particle number to 
demonstrate that our hydro-kinetic approach agrees with the main features of particle emission in 
the midrapidity region of relativistic heavy ion collisions, known from results of many transport 
code simulations. We also compare the particle momentum spectra from HKM with ones obtained 
by means of the Cooper-Frye prescription with standard isothermal freeze-out hypersurface. 

We assume the following Bjorken-type initial conditions at tj = 1 fm/c for HKM calculations: 
initial longitudinal flow vl = z/t without transverse collective expansion, ideal Boltzmann gas 
of 7r~ in chemical (local) equilibrium as the initial particle distribution, boost-invariance of the 
system in the longitudinal direction and cylindrical symmetry with Woods-Saxon initial energy 
density profile in the transverse plane, 

e(r.,rr) = j-^^ , (47) 

exp (r^^) + 1 

where Rt = 7.3 fm, 5 = 0.67 fm are associated with the density profile of the Au nucleus, and 
the maximal energy density in the center of the system, eo = 0.4 GeV/fm^, corresponding to the 
temperature T ~ 320 MeV of chemically equilibrated pions. To find the momentum spectra in 
the HKM, the relaxation time Trei = Pot*^i/p^u^ needs to be specified. For this aim we utilize 
the expression given by Eq. (|19p that represents the rate of binary collisions for one component 
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Boltzmann gas. As for the cross section, we carry out calculations for two distinct (artificial) values: 
(T = 40 nib and a = 400 mb. We performed ideal hydro-calculations till Tf = 30 fm/c when the 
system becomes very rarefied, therefore the interactions are nearly ceased and momentum spectra 
are almost frozen. The values of u^{x), T{x) and ^{x) near midrapidity, 77 ps 0, are used then to 
calculate the emission function S{t',r',p), Eq. (f28]) . and, utilizing Eqs. ([231) and ([211) to evaluate 
the momentum spectra at hypersurface r = Tf . Note that the transformation of Eqs. ([23p and ([24p 
for evolution parameter r instead of t with corresponding substitution for the initial conditions can 
be easily done in a straightforward manner. 

The results for the pion emission density integrated over the transverse momenta px = {pt cos (j), pT sin < 
at zero longitudinal momentum, pi = 0, 

{S{x))p^^^^ j S{x,p)(fpT, (48) 

as a function of transverse radius rx and Bjorken proper time r are shown in Figs. [1] and [2] for 
(7 = 40 mb and a = 400 mb, respectively. The integrated emission function is multiplied by a 
"geometric" factor r, because of kinematics in hyperbolic coordinates: r • S{t, rx) is associated 
with the probability for particles with any transverse momentum to be emitted in the central unit 
of rapidity within the proper time interval [r^T + dr] and transverse radius [vT^rx + drx]- As 
one can see in these figures, the maximum of emission is more pronounced and occupies narrower 
space-time area for larger cross sections than for smaller ones, which is certainly a result of higher 
medium opacity for higher rate of collisions. Therefore the freeze-out is rather gradual process 
for the low {a = 40 mb) value of the cross section in the expanding gas, and one could hardly 
expect the validity of the Cooper- Frye prescription there. Figures [3] and JH show the dependence of 
the emission function S on the angle between the position and momentum vectors at the particle 
emission point in the transverse plane. 

To check the reliability of the conventional CFp applied at the isotherm T = Tf ^ m-,^ versus the 
hydro-kinetic picture of the continuous particle emission, we compared the transverse momentum 
spectrum of emitted pions in HKM with the one calculated in CFp. We performed an integration 
of the emission function over a space-time four-volume till Tf = 30 fm/c to obtain the momentum 
spectrum of pions in the HKM. The emission of particles already free at the initial moment tj was 
also accounted for in such a calculation. On the other hand, we utilized CFp for an isothermal 
hypersurface with Tf = 120 MeV and locally equilibrated distribution function on it. The value of 
the freeze-out temperature, Tf = 120 MeV, is near temperatures corresponding to maximum of the 
emission function integrated over transverse momenta for an expanding gas of pions with a = 400 
mb. Such a hypersurface, as well as the velocity and chemical potential distributions on it were 



extracted from a pure hydrodynamic calculation. The results for Tn^-spectra {niT = \/pt + "^^) 
are presented in Fig. [5l One can see the increase of the effective temperature (inverse slope) in 
the case of small cross section {a = 40 mb) compared to pure hydrodynamic calculations with 
Cooper-Frye prescription. This happens because, if the collision rate is small, particles can escape 
easily from the early stages of the evolution when T > Tf (see Fig. [3]). While the transverse 
collective velocity for particles that escape at higher temperatures T > Tf is smaller than the one 
for particles which suffer freeze out according to CFp at T = Tf, the resulted gain in collective 
velocity does not lead, typically, to an increase of the effective temperature calculated according 
to ideal hydro-equations because of the energy transfer from transverse to longitudinal degrees of 
freedom. Then the utilization of CFp with Tf ~ m^ leads to a noticeable decrease of the inverse 
slope compared to the hydro-kinetic result (see Fig. |5|). Apparently, the high collision rate that 
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occurs for a = 400 mb prevents an early escaping and, consequently, leads to a gradual freeze-out in 
the rather dilute medium at low temperatures T < Tf. The interactions at this stage do not change 
the momentum spectrum essentially justifying, thereby, the utilization of CFp for the Tf = 120 
MeV isotherm. 

Now, let us consider a more sophisticated hydro-kinetic model accounting for some realistic 
features of the evolution of fireballs created in ultrarelativistic heavy ion collisions. We focus on 
the midrapidity hadronic emission at RHIC energies. Note that resonance decay contributions 
to pion spectra are not taken into account in the present version of HKM. Therefore, we do not 
compare our results with experimental data. Our aim here is to study the influence of different 
types of initial conditions and equations of state on hadronic emission processes. Special attention 
is paid to the lifetime of the system that undergoes a phase transition, because, as is well known, 
the long evolution time is one of the main shortcomings of hydrodynamic and kinetic models aiming 
to describe the quark-gluon plasma (QGP) to hadron gas (HG) transition in relativistic heavy ion 
collisions: it leads to an overestimate of the the corresponding lifetime scale found in HBT analysis 

m- 

First of all, to solve the relativistic hydrodynamic equations, an equation of state needs to be 
specified. Because the thermal model analysis of the particle number ratios at RHIC demonstrates 
almost zero baryon chemical potential [3T], we utilize an EoS with zero net baryon density. Ac- 
cording to lattice QCD calculations, the deconfinement transition at a vanishing baryon chemical 
potential is a rapid crossover rather than a first order phase transition with singularities in the 
bulk thermodynamic observables [32j. Therefore for the calculations presented here, we use for 
high temperatures a realistic EoS p3j adjusted to the QCD lattice data with a crossover transi- 
tion at about Tc ~ 175 MeV and matched with an ideal chemically equilibrated multicomponent 
hadron resonance gas at T = T^. For temperatures in the interval Teh < T < Tc, Teh = 160 MeV, 
we utilize an EoS [33^j of an ideal chemically equilibrated multicomponent hadron resonance gas, 
and for T < Teh we obtain and utilize an EoS of a multicomponent ^34j hadron resonance gas with 
a chemical composition frozen at T = Teh- The resonance mass spectrum extends over all mesons 
and baryons with masses below 2.6 GeV [M]; corresponding electronic tables where taken from the 
FASTMC event generator [35]. The value of the chemical freeze-out temperature, Teh = 160 MeV, 
is chosen because it is near chemical decoupling temperatures |31] extracted from the measured 
hadron abundance ratios at RHICo 

Then, to calculate the equations of state p = p{T, fii, ...,/i„), e = e(T, /^i, ...,/U„) of the chemically 
frozen ideal hadron resonance gas, one needs to know how the chemical potential Hi of each particle 
species and energy density depend on temperature when the hadron gas, being in an initially 
known state (all Teh,eeh,^ch ^^^ known), expands adiabatically. It can be done in the following 
way. Firstly, note that the concentration of any particle species "i" , 



ni{T,fii 



n 



Y,nu (49) 



(50) 



n(T,//i,...,/i„)- 

is fixed by its initial value at Teh and, so, 

ni{T,fii) _ 
nj{T,fij) 

Note to avoid the misunderstanding that the decay of resonances during the hydrodynamic expansion of hadronic 
matter agrees with the chemical freeze-out concept but could influence an EoS and particle spectra, we neglect here 
this effect just for simplicity reasons. 
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does not change with temperature for T < Teh- Here rii is the number density of particle species 
"i" and, in Boltzniann approximation, 

ni = — ^-"2 — Tmfexp{fXi/T)K2{mi/T), (51) 

where Kn{u) = ^ J_^ dz exp[—u cosh z + nz] , with Reu > 0, is the modified Bessel function of 
order n (n = 0, 1, ...). The energy density and pressure of the chemicahy frozen mixture of ideal 
Boltzmann gases are 



where 



n{T,fii,...,fin)'^Kiei{T), p{T, fii, ...,fin)T, (52) 



_ ei Ki{mi/T) 
Ci = — = rrii——- -— - + 3T. (53) 



Then, taking into account the thermodynamic identity 






e + p = sT + 2^ Him (54) 

and, accounting for the constancy of entropy density to particle number density ratio, s/n, during 
an ideal hydro-evolution, we divide the above identity by {nT) and get 

Equations ([5^ and ([50]) . ([^T|) define fii{T) for all particle species "i" and, thereby, complete the 
thermodynamic trajectory (for T < Teh) of the system and all their components: e{T) , p{T) , miT) 
in isentropic expansion starting from Teh- Of course, the equations of state of the chemically frozen 
hadron resonance gas (j52p can be reduced now to p(e) on this particular thermodynamic trajectory. 

For comparison, we also perform ideal hydro-calculations with the EoS taken from Ref. [36], 
where a strong first order phase transition is assumed. The corresponding EoS are presented in 
Figs. El El 

To calculate the pionic emission function and momentum spectra by means of HKM, the rate 
of collisions needs to be specified. We use the expression for the collision rate (in the co-moving 
system) that accounts for binary collisions of pions with any "i" hadronic species. 



' -y fd'kj^e.p(-^'''-/f^^)a.is)'^^'-^'^-'^'^"^^'-^'^^'^'^"\ (56) 



rltix,p) ^J (2vr)3 "V T{x) J '' ' 2E,Eu 



Here gi = (2jj + 1), Ep = y^p^ + m?^ E^^i = yk^ + fnf, s = {p + k)"^ is the squared cm. energy of 
the pair, and cri{s) is the total cross section in the corresponding binary collision. For the latter, 
we utilize a Breit-Wigner resonance formula with -y/s-dependent parametrization of decay widths 
as in Ref. [37]. All relevant resonance states from fM], used in the FASTMC event generator |35j 
- 359 different species - are taken into account for the calculation of cjj(s). Suppression of bulk 
pionic emission from very high energy density stage of matter evolution where quark-gluon degrees 
of freedom are dominated is assured in our approach by abrupt increase of the collision rate (j56p 
that is conditioned by a drastic increase of hadronic density at T > Tc in Eq. ()56p . Note, to avoid 
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a misunderstanding, we utilize the relaxation rate (j56p to calculate the pionic emission for T > T^ 
for the sake of simplicity only. In a more advanced approach the gradual disappearance of pions 
and other hadronic degrees of freedom during the crossover transition as well as the scattering of 
pions with quarks, gluons, etc., should be taken into account. Here we just expect that it will lead 
to the same effect of drastic suppression of the pionic emission from the region occupied by QGP. 

As for the initial conditions, we first perform calculations with a Woods-Saxon initial energy 
density profile (I47p in the transverse directions at tj = 1 fm/c. The values of parameters we choose 
for illustrative calculations are eo = 6 GeV/fm^, Rt = 7.3 fm, and 5 = 0.67 fm. It results in the 
initial temperature at the central "plateau" T = 247 MeV. These initial conditions at midrapidity 
are similar to the ones used in Ref. [38j, where the transverse energy density distribution at r, = 1 
fm/c was parametrized by a flat region with Gaussian smearing near the edge. The hydro-evolution 
starts at Tj = 1 fm/c with initial Bjorken flow without transverse velocity {vt = 0), assuming a 
complete chemical equilibrium at T > Teh and a chemically frozen evolution below Teh- 

The px-integrated pion emission function at midrapidity is shown in Fig. [9] for an EoS with 
a crossover transition and in Fig. [10] for an EoS with strong first order phase transition taken 
from Ref. [36]. One can see that the utilization of the EoS with a crossover transition leads to a 
significant decrease of the system lifetime because of a more effective acceleration of the system [39] 
in expansion and, therefore, is more reliable in view of RHIC data indicating a rather short lifetime 
of the system [30] [^ Results for a non-px-integrated emission function are presented in Fig. [TTl 
They demonstrate the tendency of high-pj- particles to be emitted early from the periphery of the 
system, whereas low-px particles are mostly emitted at the late stage of evolution from the center 
of the system when the system becomes fairly rarefied because of the expansion. It corresponds 
to results of Ref. [13], where the same conclusion was drawn based on partonic cascade model 
calculations. 

While the system lifetime is reduced if the EoS with a crossover transition is utilized, it could still 
be too high for reproducing measured values of the HBT radii at RHIC [30]. To study the influence 
of the initial conditions on the lifetime of the system, we perform hydro-kinetic calculations also 
for a Gaussian initial energy density profile, instead of a Woods-Saxon one: 



€(Ti, tt) = eo exp ( --^ j , (57) 

where Rt = 7.3 fm and the normalization constant is eo = 6.0 GeV/fm^, so that the total energy in 
any rapidity slice remains the same as for the Woods-Saxon initial density profile. The transverse 
momentum integrated pion emission function at midrapidity is shown in Fig. [T2l One can readily 
note a considerable decrease of the "averaged" emission time due to a faster transverse velocity 
development because of the initial gradients of density in the whole transverse region and, therefore, 
fast cooling. Also we checked the validity of the Cooper- Frye prescription for spectra at two assumed 
freeze-out hypersurfaces corresponding to the isotherms T = 75 MeV and T = 160 MeV (see Fig. 
13]) . here the value of the freeze-out temperature, Tf = 75 MeV, is chosen because it is near 
temperatures corresponding to a maximum of the emission function integrated over transverse 
momenta. We find that the effective temperature of the spectrum approximately coincides with 
the one calculated according to CFp at hypersurfaces T = 75 MeV for px < 0.5 GeV and T = 160 
MeV for pt ^ 1 GeV respectively (see Fig. [T^ . The spectrum calculated with HKM is concave, as 



^The same conclusion was reached recently in Ref. 40 based on the analysis of effects of different forms of the 
sound velocity function on the hydrodynamic evolution of matter created in ultrarelativistic heavy ion colhsions. 
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shown in Ref . [13] , and demonstrates that the utihzation of a simple p- independent isothermal CFp 
for modeling freeze-out for such systems could seriously underestimate the effective temperature of 
the momentum spectra. 

Next, in addition to the Gaussian energy-density profile, we consider also a (pre-equilibrium) 
nonzero initial flow in the transverse direction. We set 

?;T = tanh(0.3--^), (58) 

Rt 

and recalculate the normalization of energy density profile, which is eq = 5.32 GeV/fm'^, so that 
the total energy of any rapidity slab of matter in midrapidity remains the same. The results in 
the shape of the emission function are shown in Fig. [151 One can see that the inclusion of the 
initial transverse flow leads, as expected, to an even faster expansion of matter and, thereby, to 
more reduction of the mean emission time. 

5 Conclusions 

We have developed the formalism of the hydro-kinetic model [15] intended for a detailed study 
of the space-time picture of hadronic emission from rapidly expanding fireballs in A+A collisions 
and, in this way, for an evaluation of the observed particle spectra and correlations. The approach 
developed is consistent with Boltzmann equations and conservation laws, and accounts also for the 
opacity effects. 

Our analysis and numerical calculations show that the widely used phenomenological Landau/ 
Cooper-Frye prescription for the calculation of pion (or other particles) spectra is too rough if 
the freeze-out hypersurface is considered as common for all momenta of pions. The Cooper-Frye 
formula, however, could be applied in a generalized form, accounting for the direct momentum 
dependence of the freeze-out hypersurface (t{p)^ corresponding to the maximum of the emission 
function S{t{r,p),r,p) at fixed momentum p in an appropriate region of r. If such a hypersurface 
cr(p) is found, the conditions of applicability of the Cooper-Frye formula for given p is that the 
width of the maximum, which in the simple cases - e.g., for one component system or at domination 
of elastic scatterings - is just the relaxation time (inverse of collision rate), should be smaller than 
the corresponding temporal homogeneity length of the distribution function. 

The first approximation within this method is done for fireballs undergoing 3D azimuthal sym- 
metric Bjorken-type relativistic expansion. For realistic cross sections, the standard CFp, related 
to the fixed isotherm, underestimates the effective temperature of the observed spectra since, dur- 
ing an ideal hydrodynamic expansion, the effective temperature of hadronic transverse spectra is 
typically decreasing, but particles with relatively high momenta can escape from the system at an 
early hot stage. The latter process is accounted for by the proposed method. 

The recent HBT data at RHIC energies restrict a possible maximal value of the system lifetime, 
that is a serious problem for hydro-kinetic "hybrid" models ^HJ of matter evolution in relativistic 
heavy ion collisions. We studied the effects of the type/order of phase transition implemented in 
the EoS, as well as of different kinds (Woods-Saxon and Gaussian) of an initial transverse profile, 
on the duration of pionic emission. We found that a realistic EoS, motivated by lattice QCD, 
produces a faster transverse expansion, reducing therefore the lifetime of the system as compared 
with the results based on the EoS with first order phase transition. The most serious reduction 
of the lifetime is observed for the initial Gaussian density profile, because for this profile there 
are initial pressure gradients over the whole transverse slice of the system already at the initial 
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moment. It is worth noting that the Gaussian- type energy density profile naturally appears [H] in 
the Color Glass Condensate representation of the initial state of colliding nuclei (for review see, e.g., 
Ref. [42]) in high energy collisions. Also the inclusion of the initial nonzero transverse collective 
velocity leads to faster expansion of the matter and, thereby, provides an additional decrease of 
the emission time. We studied the effects of inclusion of such flow in the initial conditions for 
hydro-kinetic model calculations and found that a moderate initial flow developed at the very early 
prethermal stage of the evolution of finite partonic systems into vacuum [l3] results in a noticeable 
decrease of the system lifetime. It could indicate a way for resolving the HBT puzzle in the hydro- 
kinetic approach. Further developments of the hydro-kinetic approach and description of the data 
will be the subject of a follow-up work. 
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Figure 1: Space-time dependence of the emission function integrated over transverse momenta 
at PL = 0, for an expanding gas of pions with cross section 40 mb, initially with longitudinally 
boost-invariant flow and Woods-Saxon energy density profile in the transverse plane. 
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Figure 2: Same as in Fig. [T]but for a cross section equal to 400 mb. 
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Figure 3: Emission function for an expanding gas of pions with cross section 40 mb at r = 2 
fm/c as a function of the angle between the position and transverse momentum vectors of the 
escaping particle at the distance ry = 6 fm from the axis for different absolute values of transverse 
momentum. 
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Figure 4: Same as in Fig. [3] but for cross section 400 mb, r = 4 fm/c and rx 



fm. 
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Figure 5: Transverse spectra of pions escaped until r = 30 fm/c from an expanding fireball calcu- 
lated in HKM, with cross sections 40 mb and 400 mb, compared to the spectrum according to CFp 
applied to the local equilibrium distribution at hypersurface Tf = 120 MeV. 
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Figure 6: Pressure versus energy density for Laine and Schroder EoS with crossover transition; 
for Hirano and Tsuda EoS with strong first order phase transition; and for chemically frozen ideal 
hadron resonance gas. 
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Figure 7: Temperature versus energy density for Laine and Schroder EoS with crossover transition; 
for Hirano and Tsuda EoS with strong first order phase transition; and for chemically frozen ideal 
hadron resonance gas. 
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Figure 8: Chemical potential of pions (in the Boltzmann approximation) versus energy density for 
Laine and Schroder EoS with crossover transition; for Hirano and Tsuda EoS with strong first order 
phase transition; and for chemically frozen ideal hadron resonance gas. 
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Figure 9: Space-time dependence of the pionic emission function integrated over the transverse 
momenta for an expanding fireball with crossover transition, initially with longitudinally boost- 
invariant flow and Woods-Saxon energy density profile in the transverse plane. 
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Figure 10: Same as in Fig. [9] but for a strong first order phase transition. 
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Figure 11: Space-time dependence of the pionic emission function of an expanding fireball with 
crossover transition and Woods-Saxon initial energy density profile in the transverse plane, at 
= and different values of pt: (a) pt = 0.3 GeV, (b) pr = 0.6 GeV, (c) pr = 1.2 GeV. 
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Figure 12: Space-time dependence of the pionic emission function integrated over transverse 
momenta for an expanding fireball with crossover transition, initially with longitudinally boost- 
invariant flow and Gaussian energy density profile in the transverse plane. 



28 





Figure 13: Isothermal hypersurfases of an expanding fireball with crossover transition and Gaussian 
initial energy density profile in the transverse plane: (a) Tf = 75 MeV, (b) Tj = 160 MeV. 
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Figure 14: Transverse spectrum of pions escaped until r = 30 fni/c from an expanding fireball with 
crossover transition calculated in HKM versus spectra calculated according to CFp applied to the 
local equilibrium distribution at Tj = 75 MeV and Tf = 160 MeV. 
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Figure 15: Same as in Fig. [12] but with an initial nonzero transverse flow added. 
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